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We report on the first measurement of J/tp production as a function of charged particle 
pseudo-rapidity density dNch/drj in pp collisions at ^/s = 7 TeV with the ALICE exper- 
iment at the LHC. 3/ip mesons are detected down to pT = via their decays into e"'"e~ 
pairs at mid-rapidity {\y\ < 0.9) and into jJ.^ pairs at forward rapidity (2.5 < y < 4). 
dNch/drj is measured within \ri\ < 1. We compare results in the two different J/i/i rapidity 
ranges. Preliminary PYTHIA simulations are also presented. 

1 Introduction 

The production mechanism of heavy quarkonium states (eg. J/ijj) is very complex and is not 
fully understood. Various models such as the Color Singlet, nonrelativistic QCD approach 
(NRQCD) and the Color Evaporation Model aim to explain how a heavy resonance state can 
be produced in a hard process [H [2]- This field is very active in theory development. In 
particular, describing the J/ip production cross-section and polarization is a challenge for most 
of the models [5111111], including also the last LHC data [5J [71 [HI H] . Furthermore, J ftp production 
could be accompanied by a hadronic activity (hadrons produced in a cone around J/V") and 
it was pointed out that new observables are needed to constrain models |10j . In addition, it 
was proposed in [lTl[T2l[T3] that intial state effects could modify J ftp production due to gluon 
density fluctuations and a special transverse structure of the nucleon. 

To look at exclusive final states and not only inclusive ones, it is needed to have a full 
description of hard processes in a complete event. A description of the interplay between the 
hard and the soft components of the event as well as of color fiow and energy conservation 
is mandatory. In high energy proton-proton collisions, the total event multiplicity can have 
a substantial contribution from Multi-Parton Interactions (MPI). With MPI, several parton- 
parton interactions can occur in a single pp collision. MPI are commonly used to describe the 
soft underlying event but can also contribute on the hard and semi-hard scale, this contribution 
becoming more and more relevant with increasing energy [141 1151 116j . The NA27 experiment 
performed a study that related open charm production and underlying event properties for pp 
collisions at y/s ~ 27 GeV. It was found that charged particle multiplicity distribution in events 
with and without charm production differs by 20% |17j , indicating already a different behavior 
in multiplicity distribution. 
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In addition to MPI, other non-trivial effects could have an impact on multiplicity dependence 
and hard interactions. The measured charged particle multiplicty in pp collisions at LHC 
energies reaches values higher than in peripheral Cu-Cu at RHIC at -y/sjviv = 200 GeV [15] , 
With the high energy density reached in pp collisions at LHC energies, some models predict the 
occurence of collective behaviour at LHC energies [121 [2U], and one could consider the possibility 
of a modification of J /if] yield in high multiplicity pp events due to collective phenomena [21] . 

In this article, we report the measurement of relative i/ip production (dA'j/^, / Ay) / 
<ANj/^/Ay > at mid {\y\ < 0.9) and forward (2.5 < y < A) rapidities as a function of the 
relative charged particle multiplicity (diVch/d7y)/<d7Vch/d7/ > for pp collisions at ^/s = 7 TeV 
at LHC measured by the ALICE experiment |22j . We will then discuss a preliminary comparison 
with PYTHIA 6.4 simulations. 

2 J/i/j production as a function of the relative charged 
particle multiplicity 

A detailed description of the ALICE setup can be found in [23 . The J/ip is detected in the 
di-electron channel with the central barrel {\ri\ < 0.9) and in the di-muon channel with the 
muon spectrometer (—4 < rj < — 2.5^. The description of the part of the ALICE setup used 
in this analysis can be found in [55]. 

The results shown here were obtained by analyzing pp collisions at y/s = 7 TeV collected in 
2010. A sample of 3.0 x 10^ minimum bias (MB) events and 6.75 x 10® /i-MB triggered events 
were used for 3/ip measurement in the di-electron and di-muon channels. This corresponds to 
an integrated luminosity of 4.5 nb~^ and 7.7 nb~^ respectively for the di-electron and di-muon 
channels. The MB pp trigger is defined with a signal in one of the two VZERO detectors, plus 
one readout chip signal in the Silicon Pixel Detector (SPD), in concidence with proton bunches 
from both sides of the interaction region. The /i-MB trigger requires, in addition to the MB 
trigger, the detection of at least one muon with > 0.5 GeV/c in the acceptance of the muon 
arm. More details on the data sets, triggers, running conditions and relative normalization can 
be found in [22| . Events with an interaction vertex not within \zytx\ < 10 cm are rejected. 
Pile-up events are identified by the presence of a second interaction vertex reconstructed in 
addition to the main vertex. They are rejected if the distance along the beam axis between 
the two vertices is larger than 0.8 cm, and if the second vertex has at least three associated 
tracklets. A tracklet is defined as any combinations of two hits in the SPD layers, one hit in 
the inner layer and one in the outer. 

The charged particle density diVch/d?/ is estimated using the number of tracklets iVtrk re- 
constructed from hits in the SPD. Using simulated events, it was verified that A^trk is pro- 
portonal to dNch/drj. Fig. [I] shows the distribution of the relative charged particle den- 
sity {dNch/dr])/<dNch/dri > after correction for SPD inefficiencies. <dNch/dr] > was mea- 
sured for inelastic pp collisions with at least one charged particle in jTyj < 1 and is equal to 
6.01 ± 0.01(stat.);to i2(syst.) [24]. Vertical dashed lines show the limit of the 5 bins in multi- 
plicity used in this analysis. 



^In the official ALICE reference frame the muon spectrometer is located at negative z positions and thus 
negative (pseudo-)rapidities. Since pp collisions are symmetric relative to y = 0, we have dropped the minus 
sign when rapidities are quoted. 
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Figure 1: Distribution of the relative charged particle density (diVch/d77)/<dA^ch/d77 > at mid- 
rapidity (|?7| < 1.0) after correction for SPD inefficiencies. Vertical lines indicate the boundaries 
of multiplicity bins considered in this article. 



For 3/^ measurement in the di-electron channel, tracks are selected by requiring transverse 
momentum > 1 GeV/c and a pseudo-rapidity cut of \ri\ < 0.9. Particle identification is 
performed by measuring the specific energy deposit dE/dx in the Time Projection Chamber 
(TPC) of the central barrel. The invariant mass distributions of e"*"e~ pairs are measured in 
intervals of the charged particle mulitiplicity as measured via the SPD tracklets. Examples of 
such mass distributions are shown, for the lowest and highest multiplicity intervals in the two 
left panels of Fig. [2] The track rotation method (green squares in Fig. [2j left panel) is used 
to describe the combinatorial background in each multiplicity intervals as well as the like sign 
distributions (open blue circles in Fig. [2j left panel). 

For J/?/' analysis in the di-muon channel, muon candidates are selected by requiring that at 
least one of the two muon candidates matches a trigger track reconstructed from at least three 
hits in the trigger chambers of the muon spectrometer. To remove muons produced at small 
angles that have crossed a significant fraction of the beam shield, a cut is applied on the radial 
coordinate of the track at the end of the front absorber (i?abs >17.5 cm). To reject events very 
close to the edge of the muon spectrometer acceptance a cut is applied on the rapidity of the 
pair (2.5 < ?/ < 4). To obtain the number of J/?/; in each multiplicty interval, a fit is used 
on the corresponding di-muon invariant mass distribution in the range 2 < Mi„v < 5 GeV/c^. 
The line shapes of the i/ip and V'l^S) are parametrised using a Crystal Ball function, while 
the underlying continuum is fitted with the sum of two exponential functions. Details on the 
quality of the fit results can be found in [6] . The two right panels of Fig. [2] show the measured 
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Figure 2: Opposite sign invariant mass spectra of the selected e~^e~ pairs [(a),(c)] and n~^iJ,~ 
pairs [(b), (d)] (full symbols) for the lowest [(a),(b)] and highest [(c), (d)] multiplicity bins. The 
number of events quoted in the figures refers to the corresponding minimum bias events. 



di-muon invariant mass distributions together with the results of the fit procedure for the lowest 
and highest multiplicity intervals. 

Fig. [3] presents the ratio of the J/ip yield in a given multiplicty interval relative to the min- 
imum bias yield. Corrections regarding geometrical acceptance, reconstruction efficiency and 
their systematics cancel out in the ratio {dN,y^/dy)/<dNjf^,/dy >. It was checked by Monte 
Carlo simulations that these corrections do not depend on diVch/d?/, in the range under consid- 
eration {dNch/drj < 32.9). The number of events used for the normalization of <dN,y^,/dy > is 
corrected for the fraction of inelastic events not seen by the trigger condition. After applying ac- 
ceptance and efficiency corrections these values correspond to the values than can be extracted 
from data published in [S]: <dNj/^/dy >= (8.2 ± 0.8(stat.) ± 1.2(syst.)) x 10^ for J/?/; ^ e+e" 
in \y\ < 0.9, and <diVj/^/d?; >= (5.8 ± 0.2(stat.) ± 0.6(syst.)) x 10^ for in 
2.5 < y < 4. 

For the di-electron analysis, the uncertainty due to background subtraction was obtained 
as the absolute differences using the like-sign and the track rotation methods. It is found to 
be between 2% and 12% for the different multiplicity intervals. For the di-muon analysis this 
is evaluated by varying the functional form of the background description (polynomial instead 
of exponential). It depends on the signal to background ratio and varies between 3% and 4%. 
For the muon measurement an additional systematic uncertainty comes from pile-up events 
and is estimated to be 6% in the first multiplicity interval and 3% in the others. To account 
for the possible changes of the pt spectrum with event multiplicity, an additional systematic 
uncertainty is determined by varying the < > of the 3/ip spectrum that is used as input to 
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Figure 3: dNj/^/dy as a function of the charged particle multipHcity density at mid-rapidity 
dNch/dr]. Both values are normalized by the corresponding value for minimum bias pp collisions. 
Two measurements are shown: in red circles at mid-rapidity {J/ip — > e+e~, < 0.9) and in 
blue squares at forward rapidities {J/ip M^/^ ^ 2.5 < y < 4). Statistical uncertainties of 
the J/ip yields are represented by error bars, while boxes reflect the quadratic sum of the 
point-to-point systematic uncertainties. 



the determination of the Monte Carlo corrections between 2.6 and 3.2 GeV/c. A systematic 
uncertainty of 1.5% (3.5%) is found for the di-electron (di-muon) analysis. The total systematic 
error on {dNj/^/dy)/<dNj/^/dy > is given by the quadratic sum of the different contributions 
and amounts to 2.5 — 12 % depending on the multiplicity interval for the di-electron result. In 
the case of the di-muon analysis, it varies between 8% in the first and 6% in the last multiplicity 
interval. An additional global uncertainty of 1.5% on the normalization of <dNj/^/dy > is 
introduced by the correction of the trigger inefhciency for all inelastic collisions. More detailed 
explanations on the estimation of systematic uncertainty estimated in this analysis can be found 
in 

In Fig. [3] an approximately similar linear increase of the relative J/ip yield (dA^j /dy) / 
<dNj/^/d?/ > with {dNc\Jdrii)/<dNch/d-q > is observed in both rapidity ranges. The enhance- 
ment relative to minimum bias J/ip yield is a factor of approximately 5 at 2.5 < y < 4 (8 at 
\y\ < 0.9) for 4 times the minimum bias charged particle multiplicity. A possible explanation for 
the observed correlation could be that 3/ip is always accompanied by a strong hadronic activity 
biasing high multiplicity events. Such a mechanism could imply particular spatial distributions 
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Figure 4: Relative i/i) yield (dAfj/^/dy)/<dA'^,j/^/dy > as a function of the relative 
charged particle multiplicity density around mid-rapidty {dNi^^/dri)/ KdN^h/'^Tl > calculated 
with PYTHIA 6.4.25 in the PERUGIA 2011 tune MSTP(5)=350. Results are shown at for- 
ward (2.5 < y < A) and central [\y\ < 0.9) rapidities for J/ip produced in hard scatterings via 
the NRQCD framework (MSEL=63). 



and J/i/i-hadron correlations could clarify the situation. Another possible mechanism would be 
initial density fluctuations accompanied by a specific structure of the nucleon. This mechanism 
seems to explain a factor 4 to 5 for charged multiplicity in J/^ events 5 times larger than in 
minimum bias events |13j . 



3 First PYTHIA 6.4 comparison 

To compare our results to the predictions of a model, the model has to be able to reproduce all 
aspects of an event : the hard, the soft and semi-hard part of the event in a consistent framework. 
Event generators seem well suited to fulfill such a requirement. The considered event generator 
should also include heavy quarks (c and b) and heavy resonances such as J/V" considering the 
correct masses of heavy quarks, energy conservation and color flow. Few models are left and 
none of them was built considering this new observable of 3/ijj yield versus multiplicity. The 
first one we can think of is PYTHIA 6.4 [3S] which is a pp event generator commonly used at 
the LHC. This is not the best model adressing quarkonium production, but it is extensively 
used, tuned and debugged. In this sense, this study is a first attempt of model-comparison for 
J/ip yield versus multiplicity. 

Detailed explanations on the physics model of PYTHIA 6.4 can be found in [5S]. A pp 
event is composed by a hard 2 — >■ 2 process. The two partons from the incoming protons 
can evolve through initial state radiation processes (ISR) before the hard subprocess. The 
two produced partons can also evolve with final state radiation processes (FSR). To produce 
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J/'0, several models are available, the color singlet (ISUB=86 with gluon fusion) where a J/f/' is 
produced with a gluon to ensure momentum conservation, and NRQCD (ISUB=421-439) where 
a pre-resonant state is produced. This pre-resonant state will then decay into a J/?/' emiting 
a soft gluon. In addition to the first hard interaction, a MPI scenario with varying impact 
parameter allows other incoming partons to undergo hard and semi-hard interactions. The first 
one is the hardest one. The following ones are ordered in hardness. In the last model of MPI 
implemented in PYTHIA 6.4, other interactions also evolve through ISR and FSR. In PYTHIA 
6.4, quarkonia production in hard process is only available in the first hard interaction and not 
in the following ones of the MPI scenario. This was corrected in PYTHIA 8. In addition to 
the hard production, J/V' can be produced from B decays. This source can be easily turned off, 
imposing B hadrons to be stable. J/'f/; can also be produced by the fragmentation of a cluster 
made of a cc pair produced by the branching of one or two gluons into cc pairs. Those gluons 
come from ISR and FSR. At the end of the procedure all produced partons are connected with 
strings that will fragment into hadrons via the LUND string model i25) . 

In this study the version PYTHIA 6.4.25 is used with Perugia 2011 (P2011) tune, a recent 
tuning of PYTHIA including the fit of LHC pp data at 7 TeV. It is well suited to describe 
multiplicity distributions. We consider in this article only MSTP(5)=350 which is the major 
P2011 tune [26 . B hadrons are forced to be stable. Two different data samples are generated. 
The first one is minimum bias (MSEL=1) for the denominator of our observable. The second 
sample used in the denominator contains J/V'- It is generated with MSEL=63 : J/V' are 
generated in hard processes with the NRQCD framework implemented in PYTHIA [3S]. This 
second sample also contains J/?/; produced by cluster fragmentation. The identification of J/?/) 
from hard processes and ^ from cluster fragmentation can be done with the flavour code of 
the mother, allowing the separation of the two contributions. Finally, the same kinematical cuts 
than in the ALICE experimental data sets are used to identify J/i/i and charged particles. As 
in the data analysis the samples are normalised by the number of events and the mean number 
of charged particles. The obtained results are shown in Fig. |4]for PERUGIA 2011 tune with 
MSTP(5)=350. Results for J/-0 produced in hard scattering only are presented for the same 
rapidity regions as in the data analysis. 

Concerning cluster contribution, further studies indicate that this contribution is strongly 
enhanced by the color reconnection mechanism. Clusters come from gluons originating from ISR 
and FSR. Therefore, this contribution should scale with multiplicity, because high multiplicity 
events correspond to events with a high partonic activity (MPI, ISR, FSR) and so an enhanced 
probability to form a cluster. Nevertheless this contribution in PYTHIA 6.4.25 does not seem to 
be under control due to the way color reconnection is implemented |26| and warnings concerning 
the physics model of the cluster implementation can be found in the PYTHIA manual pSj . 
Thus, this contribution is excluded in results of Fig. [4j 

The i /i}) production from hard processes only, decreases as a function of relative multiplicity, 
which is in disagreement with the trend observed in the data (Fig. [3]) . A nai've interpretation 
would be that if the first hard process is independent from the other aspects of the event (ISR, 
FSR, MPI, underlying event), the i/ij) yields should be flat as a function of multiplicity. It is not 
observed with PYTHIA 6.4.25, and further studies with other models are needed to understand 
this behaviour. E.g. a study with PYTHIA 8 might provide further insight. In PYTHIA 8, 
MPI could play a role in J/?/; production with the possibility of quarkonium production in all 
the parton-parton interactions of the MPI scenario. Other event generators, such as CASCADE 
[27], should be tested as well. 
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4 Conclusions 



In summary, we have presented the first measurement of the J/'0 yield as a function of the 
charged particle multiplicity dNch/dr], performed by the ALICE collaboration [5^. 3/ip mesons 
are detected at mid-rapidity (|y| < 0.9) and forward rapidity (2.5 < y < 4), while dNch/drj 
is determined at mid-rapidity (|r/| < 1). An approximately linear increase of the J/ip yields 
with the charged particle multiplicity is observed. The increase is similar at forward and mid- 
rapidity, exhibiting an enhancement relative to minimum bias J/ip yield by a factor of about 5 
at 2.5 < 2/ < 4 (8 at lyl < 0.9) for ^ 4 times the minimum bias charged particle multiplicity. 
A first comparison with PYTHIA 6.4 simulations was shown. J/^ yields originated from the 
first hard interaction do not follow the same trend as seen in the data. Further studies are 
needed to explore the physics involved in such a new observable. From the experimental side 
the study of charged particle multiplicity dependence of T, open charm and also other hard 
observables such as jet and Drell-Yan production could bring more informations. Multiplicity 
studies for various classes are also of interest. One could also propose an underlying event 
study replacing the leading jet hy a J/ip and J/^-hadrons correlation studies. From the event 
generators side, new studies are needed with PYTHIA 6.4, considering new tunes and other 
sets parameter. Other models should also be tested such as PYTHIA 8 and CASCADE. 
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